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Abstract 

The quantum cosmology of two-dimensional dilaton-gravity models is in- 
vestigated. A class of models is mapped onto the constrained oscillator-ghost- 
oscillator model. A number of exact and approximate solutions to the corre- 
sponding Wheeler-DeWitt equation are presented. A wider class of minisuper- 
space models that can be solved in this fashion is identified. Supersymmetric 
extensions to the induced gravity theory and the bosonic string theory are then 
considered and closed-form solutions to the associated quantum constraints 
are derived. The possibility of applying the third-quantization procedure to 
two-dimensional dilaton-gravity is briefly discussed. 
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1 Introduction 



It is widely thought that quantum gravitational effects become important on scales 
of the order of the Planck length. It follows, therefore, that insight into the nature of 
quantum gravity might be gained by considering the very early Universe and the end- 
point of black hole evaporation. However, there remain many unresolved technical 
and conceptual difficulties with (3 + l)-dimensional quantum gravity and to make 
progress one must consider simplified toy models. 

Recently there has been considerable interest in two-dimensional theories of grav- 
ity. The Einstein action is a topological invariant in two dimensions and must there- 
fore be modified if the theory is to be non-trivial. The simplest extension is to 
include a non-minimally coupled, self- interacting scalar 'dilaton' field. These models 
are closey related to string theory in non-critical dimensions ]IJ and may provide a 
solvable framework in which some of the questions raised in quantum gravity can be 
studied. Indeed, black hole evaporation has been extensively investigated following 
the introduction of the Callan, Giddings, Harvey and Strominger (CGHS) model @. 

Quantization of two-dimensional models has been performed by employing a num- 
ber of techniques such as the BRST and Dirac operator methods ]3], |J]. However, it 
is possible that lower-dimensional gravity may also shed light on some of the issues 
raised in quantum cosmology, such as the problem of extracting physical predictions 
from the wavefunction of the Universe. In this paper we investigate the quantum 
cosmology of a generalized class of two-dimensional models. 

A number of approaches to quantum cosmology may be taken. In the path integral 
formalism, for example, the wavefunction of the Universe is expressed as a path 
integral over a certain class of metrics, matter distributions and manifolds |5], |J. 
This formalism has recently been investigated within the context of two-dimensional 
quantum gravity by Ishikawa 0. 

Alternatively, one may follow the Dirac quantization procedure ||. The wavefunc- 
tion of the Universe is annihilated by the Hamiltonian operator and, in principle, its 
functional form can be determined by solving the zero-energy Schrodinger equation. 
This equation decomposes into the Wheeler-DeWitt and momentum constraint equa- 
tions that describe, respectively, the invariance of the theory under reparametrizations 
of time and spatial diffeomorphisms [[|. 

In general, the Wheeler-DeWitt equation is a functional differential equation and 
is very difficult to solve. Solutions can be found in four- dimensional theories by 
invoking the 'minisuperspace' approximation and freezing out the inhomogeneous 
modes. It is not yet established that such an approach can lead to meaningful results 
in higher dimensions, although arguments have been developed to suggest that it may 
be relevant [|l(|. On the other hand, the existence of a Killing vector in a wide class 
of spatially closed, two-dimensional cosmologies implies that all classical solutions to 
the field equations are spatially homogeneous ||. It can therefore be argued that the 
minisuperspace approach is exact for these models 0, IT21 . 
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Adi and Solomon have adopted a geometrical approach and found a new solution 
of the Wheeler-DeWitt equation JE3]. Navarro-Salas et al. |Tl| have quantized the 



induced gravity theory |14[] via the covariant phase-space and reduced ADM phace- 
space methods. They also derived and solved the Wheeler-DeWitt equation for this 
model. Henneaux, on the other hand, performed the quantization of this theory in the 
functional Schrodinger representation by first solving the supermomentum constraint 
at the classical level . This technique was subsequently generalized to other models 

E3- 

In this paper we follow the approach normally employed in four-dimensional quan- 
tum cosmology. We consider spatially closed cosmologies within the generalized class 
of two-dimensional dilaton-gravity models, where the dilaton field is assumed to be 
constant on the surfaces of homogeneity. The models are defined by the functional 
form of the dilaton potential and specific models have been considered previously in 
this fashion by a number of authors [[□], [H], [TJ, [HJ . We find that the Wheeler-DeWitt 
equation is exactly solvable for a wide class of potentials. 

The paper is organized as follows. The classical dynamics of these models is 
considered in Section 2 and a subset is mapped onto the constrained oscillator-ghost- 
oscillator model. In Section 3, these models are quantized and a number of exact 
solutions to the Wheeler-DeWitt equation are presented. An interpretation of these 
solutions is discussed within the context of a specific model. Two classes of approxi- 
mate solutions are presented in Section 4. The first is a power series solution to the 
Wheeler-DeWitt equation derived by employing a modification of the Picard iteration 



scheme [2D]. The second is the class of WKB solutions derived by means of a Legendre 



transformation [pip . In Section 5, a wider class of exactly solvable two-dimensional 
minisuperspaces is identified. If appropriate conditions are satisfied, these models can 
also be mapped onto the constrained oscillator-ghost-oscillator model. The superpo- 
tential of the wavefunction must be a separable function of the minisuperspace null 
coordinates. We employ this observation to transform the Wheeler-DeWitt equation 
derived from a renormalizable dilaton-gravity model into a solvable form |fL9[| . In 



Section 6, a supersymmetric extension is considered. It is found that the superspace 
Hamiltonians derived from the induced gravity theory and the string effective action 
may each be viewed as the bosonic component of a supersymmetric Hamiltonian. 
This symmetry is preserved at the quantum level and the associated quantum con- 
straints are solved exactly in closed form for both models. We conclude in Section 7 
with a brief discussion on the possibility of applying the third quantization procedure 
to two-dimensional cosmologies. 

Units are chosen such that h = c = 1 unless otherwise stated. 



2 Two-dimensional dilaton-gravity models 
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2.1 The generalized action 



The most general action for two-dimensional dilaton-gravity that is invariant under 
local reparametrizations and does not contain third or higher-order derivatives is 



22 



S 



'27 ' 



-g [c 2 {$)R + ci($)<r0„$0„$ + U{<f>)\ , 



(2.1; 



where g^ u is the metric on the two-dimensional space-time manifold, g is its determi- 
nant, Cx($) and c 2 ($>) are functions of the dilaton scalar field R is the curvature 
scalar and U(§) is the dilaton potential. 

Specific forms for the functions {c\, c 2 , U} correspond to different two-dimensional 
models. For example, the induced gravity action is a special case of Eq. (|2.1| ) with 
Ci = 1, c 2 = 2$ and U = 4A 2 P3"[| . The constant curvature condition may be derived 
from this action and provides a suitable analogue to the Einstein field equations in 
two dimensions |1~4||. Also of interest is the spherically-symmetric four- dimensional 
Einstein-Hilbert action. This is equivalent to Eq. ( |2.1|) if c 2 = e~ 2<1> , c\ = 2c 2 and 
U = 2(1 — Ae~ 2 *), where A is the four-dimensional cosmological constant. In this 
example the dilaton field is related to the radius of the two-sphere. 

Actions of the form ( |2.1| ) also arise in string theory [|I|]. To leading order in the 
inverse string tension a', the tree-level effective action for the closed bosonic string 
in two dimensions is given by 



S 



/ d 2 x\ 

2vr 



-2$ 



i? + 4(V$) 2 + £>($) 



(2.2) 



where -D($) = c = 16/a' is proportional to the effective central charge and the 
tachyon field is assumed to be zero |24|] . The field strength H^x of the antisymmetric 
tensor field vanishes identically in two dimensions. This action is closely related to 
the gravitational sector of the CGHS model 0. 

Closed string loop corrections introduce additional terms into the beta functions 
and therefore modify the effective action (|2.2j) p5| . If field derivatives are neglected, 
only the dilaton potential -D($) is altered and, in general, the loop corrections have 
the form 



Dm 



n>0 



,2n$ 



where clq = c and n represents the number of handles inserted [|26 
other coefficients determined by the string theory. 



(2.3) 



The values of the 



McGuigan, Nappi and Yost have considered two-dimensional string theories con- 

They showed that open strings are governed by the Born- 
1 . Orientable open strings may couple 



taining gauge fields |27 
Infeld action for non-linear electrodynamics 
to SU(N) and non-orientable strings to SO(N) or Sp(N), where must be even 
for non-orientable strings |29[|. When hole and crosscap corrections []30| are included, 
it is found that the modified equations of motion may be derived from an effective 
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action of the form (Q[), where 

£>($) = c - k(N + 2??)e*, (2.4) 

and i] = — 1,0, +1 when the gauge group is SO(N), SU(N), or Sp(N), respectively. 
k is a positive-definite open string coupling constant. 

Mignemi has recently investigated action ( |2.2j ) with D = Ae _2h * for some arbitrary 
constant h, and has found black hole solutions |31|. This action reduces to the 
string effective action in the limit /i — ^ and is also conformally equivalent to a 
two-dimensional higher-order, pure gravity theory with a Lagrangian given by L = 

R h/(h-l)^ where h ^ x H 

The general action ( |2.1| ) may be simplified after suitable redefinitions of the dilaton 
and graviton fields fl2~2"| . If c\ and C2 are positive-definite functions, we may define a 
new scalar field 

E = V2J d$^/ci($) (2.5) 
and perform the conformal transformation 

9^ = ^V> fi2 = e_2p ; (2-6) 

where 

p=%-\r**(%Y (2.7) 



(2.8) 



q 2 A J \dE' J 

and q is a constant. The action (|2.1| ) transforms to 



2tt 

where 



5 = / d 2 xJ—g 



g^ = 2c 2 [$(S)], V(^) = e 2 "t/($). (2.9) 

Thus, the models are defined uniquely by the functional form of the dilaton potential. 
Unless otherwise stated, we shall view theory ( |2.8j ) as our starting action and we 
therefore drop the tildes for notational simplicity. 



2.2 Two-dimensional cosmologies 

In the canonical framework the topology of space-time is S x 9ft, where the real line 
3ft corresponds to the time dimension and the spatial section E is either a real line or 
a circle S 1 . The former case applies to two-dimensional black hole solutions, whereas 
the compact spatial topology is relevant for cosmological models. In this case the 
world-interval has the form 

ds 2 = -N 2 (t)dt 2 + a 2 (t)dx 2 , (2.10) 
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where a(t) is the radius of the spatial hypersurfaces and N(t) is the lapse function. 
The spatial sections represent surfaces of constant ip. When the line element is given 
by Eq. ( p. 10 ), the Ricci curvature satisfies y/—gR = —2dt(aK), where the extrinsic 
curvature scalar is given by K = —a/(aN) and a dot denotes differentiation with 
respect to t. It follows, therefore, that the action (|2.8|) takes the form 



S 



dt 



(qipd+~a'tp* ) - 



(2.11) 



after integration over the spatial sections. 

We now proceed to express the kinetic terms of Eq. (|2.11|) in canonical and diag- 
onal forms by introducing new variables. This will allow the classical field equations 
to be solved for appropriate choices of lapse function and leads to simple forms for 
the Hamiltonian constraint. We begin by defining the new coordinate pair 



■ft 



= V2qae^ /2q , 



= V2qe~ 



(2.12) 



The range of these variables is determined by physical considerations. The physically 
interesting region of parameter space is < {a, qip} < +oo and this corresponds to 
the range < u < +oo and < v < V%q. In terms of these variables the action 
( p. 11| ) has the form 

S = - J dt -^iiv + NaV (2.13) 

and the Hamiltonian constraint, derived by functionally differentiating the action 
with respect to the non-dynamical lapse function, is given by 



1 . . uv 



0. 



(2.14) 



The second term on the left-hand side is the 'superpotential'. Its direct dependence 
on the dynamical degrees of freedom may be eliminated by introducing the rescaled 
variables 

u 2 

and 



(2.15) 



13= / dv'v'V[ip(v')} 



-Q 



d^e-^' /q V(^') 



(2.16) 



where (3 represents a rescaled dilaton field. It follows that < a < +oo, but the range 
of values spanned by (3 is model-dependent. The Jacobian of this transformation 
vanishes if V(ip) vanishes and we therefore restrict our discussion to potentials that 
are either positive- or negative-definite. The action ( |2.13| ) now takes the form 



S = - dt 



NaV 



a(3 + NaV 



(2.17) 
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and the corresponding Hamiltonian constraint (|2.14|) becomes 

H = aV(p aP0 - 1) = 0, (2.18) 

where p a = —0/(NaV) and p@ = —a/(NaV) are the momenta conjugate to a and 
P, respectively. 

In the gauge iV -1 = aV, the field equations take the simple form a = f3 = and 
have the general solution 

a = a Q + b(t — to) 
P = P + b- 1 (t-t ), (2.19) 

where {a , p , t , b} are constants. We conclude, therefore, that the classical dynamics 
of these two-dimensional Universes is equivalent to that of a non-interacting point 
particle propagating on two-dimensional Minkowski space. The variables a and P 
may be viewed as null coordinates and the regime of Minkowski space accessible to 
the 'particle' depends on the dilaton potential. 

We now impose the additional restriction that ft remains positive- or negative- 
definite for all physically relevant values of the scale factor and dilaton field. In this 
case we may introduce a third pair of variables defined byQ 



w = yjjp - yfa, z= yJjP + y/a, (2.20) 
where 7 = P/\P\- For real values of a and P, z > \w\. In the gauge 



N=±-(a 1 P) 1/ \ (2-21) 
av 



action ( p. 17 ) transforms to 



S = - dt 



7 (z 2 - w 2 ) + - (z 2 - w 2 ) 



(2.22) 



This is the action for the constrained oscillator-ghost-oscillator system when P < 
and the model corresponds to a constrained hyperbolic system if P > 0. In this latter 
case, the field equations have the general solution 

w = Ae 1 ' 2 + Be- 1 ' 2 

z = Ce t/2 + De- t/2 , (2.23) 

where the constants of proportionality satisfy AB = CD and are chosen to ensure 
z > \w\ for all t. The trajectory of this solution is a central conic section in this sector 
of the (w, z) plane. If P < 0, however, the general solution is 

w = e£cos(t/2 + #i), z = Ecos{t/2 + e 2 ) , (2.24) 



1 This restriction on j3 ensures that the corresponding Whceler-DcWitt equation will not be of 
the elliptic type when expressed in terms of w and z. 



6 



where {E, 81,82} are arbitrary, real constants and |e| = 1. These solutions lie on the 
family of ellipses |[T7| 

w 2 + z 2 -2ewzcosd = E 2 sm 2 8, (2.25) 

where the eccentricity is determined by 8 = 9\ — 82 and the major axis lies along the 
line w = ez. 

This correspondence between the constrained oscillator-ghost-oscillator and a spe- 
cific two-dimensional dilaton-gravity model was recently observed by Onder and 
Tucker [|17] in the synchronous gauge N — 1. Their model corresponds to the choice 
ci = c, C2 = |$ and U = A + Ae c * in action (|2~T|), where {c, A, A} are constants. 
Within the context of this model, they employed such a correspondence to investi- 
gate the connection between the classical and quantum cosmologies. They identified 
appropriate linear superpositions of quantum states that highlighted the classical or- 
bits (|2.25| ) and were therefore able to conclude that a definite correlation between 
classical and quantum solutions exists in this model. This is interesting because the 
question of how a classical space-time emerges from a quantum theory of the Uni- 



verse is currently unresolved. The above analysis generalizes the results of Ref. [17| 
and shows that the correspondence between the oscillator model and two-dimensional 
cosmologies arises in a wide class of dilaton-gravity models. This suggests that two- 
dimensional theories may provide valuable insight into the problems associated with 
quantum cosmology in higher dimensions. In view of this we proceed in the next 
Section to investigate the quantum cosmological behaviour of these models. 



3 Exact quantum wavefunctions 



3.1 The Wheeler-DeWitt equation 



The cosmological models defined by Eq. (|2.17|) are quantized with the algebra 
[a,p a ]_ = i and [f3,pp\- = i. The Wheeler-DeWitt equation is the operator form 
of the Hamiltonian constraint (|2.18| ) and is realized by identifying p a with —id /da 
and pp with —id/ 8(3. The physical states, of the Universe are annihilated by this 
Hamiltonian operator. We shall not consider the ambiguities that arise in operator 
ordering, so the Wheeler-DeWitt equation has the form 



2 



dad(3 



+ 1 



^ = 0. 



This equation admits a number of exact solutions One family is given by 



-iba—i(3/b 



(3.i; 



(3.2) 



where b is an arbitrary, complex constant, j^j is bounded everywhere when Im& = 
and is bounded for Im& < if (3 < 0. 
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A natural generalization of this solution is to include a variable amplitude A(a, (3). 
Substitution of the ansatz ^ = A^/ b into Eq. implies that A satisfies 



(3.3) 



3 2 A _ t OA _ ib dA 

dad(3 b da d(3 

and one non-trivial solution to this equation is 



A = ba - 



(3.4) 



In terms of the coordinate pair fl2.2U| ) the Wheeler-DeWitt equation ( |3.1| ) trans- 
forms into O, |33], |34], p5l BO, B7I 



2 



d 2 



2 2 
W — Z 



* = 0. 



<9u> 2 <9;z 2 

This has separable solutions of the form $ = X^c^n, where 

3 -cU 2 +2 2 l/2 



' r 



a/c' 



y/cz 



(3.5) 



(3.6) 



c n are arbitrary complex coefficients, H n is the Hermite polynomial of order n and 
c = z (c = 1) for 7 = +1 (7 = —1). 

If (3 does not change sign, a third class of solution is generated by defining the 
variables 



s = - In (Aaryp) 
o 

r = - In \ — I 

6 

and the Wheeler-DeWitt equation becomes 

d 2 d 2 



(3.7) 



ds 2 dr 



, + 9 7 e 6s 



* = 0. 



(3.8) 



The wavefunction is an eigenstate of the momentum operator d/dr and has the sep- 
arable form 



e ipr Z ±ip/3 v^e 



.3s 



(3.9) 



where p is a separation constant and Z± ip / 3 represents a linear combination of ordinary 
Bessel functions of order ±ip/3. 

It should be noted that technical questions arise when quantizing with variables 



that are restricted to a finite range, as is the case in the derivation of Eqs. (|3~1 



and (^.5|). However, these issues are beyond the scope of the present work On 
the other hand, the variables ( |3.7| ) are unrestricted and Eq. ( |3.8| ) can be derived 
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from the corresponding classical action with an appropriate choice of factor ordering. 
Eqs. (|3.1|) and flO] ) may then be derived directly from this equation by a change of 
variables. 

An additional class of exact solutions may be generated by defining new variables 



a i — - a 



fi = - + y/2j/3, v = - - p^. (3.10) 



The Wheeler-DeWitt equation transforms into 

" d 2 d 2 

+ 7(/i- z/) 



^ = o (3.11; 



9/i 2 dv 2 
and has the separable solution 

$ = \c\K\{m — 7/i) + c 2 Bi(m — 7/2)] [c 3 Ai(m — 71/) + c 4 Bi(m — jv>)] (3-12) 
in terms of Airy functions, where {m,Cj} are arbitrary constants p9| . 

3.2 Exponential dilaton interactions 

It is useful to consider a specific model in order to discuss these solutions. We shall 
investigate an action of the form fl2.2|) , where the dilaton potential is given by D = 
Ae~ 2M> for some constants {h, A} |31|. This action is conformally equivalent to Eq. 
where 

# = 2e" 2 *, VW = j(^e+l*. (3.13) 



The kinetic terms in the action (|2.11|) may be diagonalized by introducing the vari- 
ables 



x = JL ( ae v-/2 9 _ e -v/2^ (314) 



and this implies that the minisuperspace metric is transformed into the Minkowski 
metric, where T is the timelike coordinate and X is the spacelike coordinate. Since 
T > only the interior region of the future light cone of the origin is covered 
by this coordinate system. We deduce, therefore, that the variables (|2.12|) represent 



the null coordinates u = T + X and v = T — X in this region. It should be noted 
that only a finite region of the interior of the future light cone is represented because 
v < fmax = V%q is bounded from above. 

For this model the (at, 0) coordinates are given by 

a = \a 2 e^ lq (3.15) 
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and 

41 + h 

/3 = -^hn», ^ = "1, (3-16) 

respectively, and when h ^ —1, the sign of /3 is uniquely determined by the sign of 
A/(l + /i). 

The parameters {a,/?} may also be viewed as null coordinates over a region of 
Minkowski space spanned by the timelike coordinate p = a + (3 and spacelike coor- 
dinate q = a — (3. Thus, if (3 > 0, it follows that p > \q\ and the analysis is again 
restricted to the interior of the future light cone of the origin. Since a, (3 G (0,oo), 
the whole of the interior is now covered. On the other hand, if (3 < 0, q > \p\, and 
the propagation of the wavefunction is restricted to the Rindler wedge of Minkowski 
space. 

The simplest interpretation of the wavefunction identifies an oscillating solution to 
the Wheeler-DeWitt equation as a Lorentzian geometry and a cosmological singularity 
is associated with an infinite number of oscillations P0|. A non-oscillating solution 



represents a classically forbidden Euclidean geometry. Let us consider the case where 
(3 < and h — 1. Eq. fl3.6|) represents the basis for a discrete spectrum of Euclidean 
solutions, where the parameter n determines the excitation level of the wavefunction 
36fl . The ground state is associated with n = and excited states with n > 0. This 



ground state is identical to solution ( |3.2| ) when b = —%. Hence, we may view the 
solution ^b=-i as the ground state of a continuous spectrum of excited states (|3.2| ) 
that are parametrized by the separation constant b with Re6 = and Im& < 0. 

Although these classes of Euclidean solution appear to correspond to classically 
forbidden behaviour, Lorentzian wavefunctions may be generated from appropriate 
linear combinations of the excited states. A more general solution to Eq. ( |3.1| ) is 
given by |33 

^= / dcM(c)e- ca -^ /c , (3.17) 



c 



where M(c) is an arbitrary function of the parameter c = ib and C represents some 
contour of integration in the complex plane. If M(c) = ^c^ ip ~ 3 ^ 3 and the contour 
of integration is over the positive half of the real axis, Eq. (|3.17 ) may be evaluated 
exactly in terms of the modified Bessel function: 

* p = lr dcciW ~ 3)/3e ~ Ca ~ lt3l/C = e lWK ip/3 (2 vW) • (3-18) 



We recognize this superposition as solution (33) with Z = K(e 3s ). This solution may 
also be generated from a linear combination of harmonic oscillator wavefunctions and, 
in general, solutions (|3.2|), (|3.6p and (|3.9p may be expressed as linear combinations 



of one other for positive and negative (3 |3| 
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The modified Bessel function has the asymptotic form K q {x) oc x~ q for \x\ <C 1 and 
q 7^ 0. Thus, the wavefunction has the form ^ p oc e ip( - r±s ^ for small spatial geometries 
(s — > — oo) and these represent plane waves in the variables (r, s). The wavefunction 
oscillates an infinite number of times when the spatial volume of the Universe vanishes 
and we identify this point as a cosmological singularity. However, the wavefunction is 
exponentially damped for e 3s > \p\/3 and this region of minisuperspace is classically 
forbidden. It is interesting to relate this solution to the classical solution in terms 
of the variables ( |3.7| ). The gravitational field equations derived from action (|2.17|) in 
the gauge N = 9e 6s / (AaVj) are given by 



f = 0, 



27 
s = —e 
47 



6s 



(3.19) 



and the general solution satisfying the Hamiltonian constraint is 



±t 



3A 



-3s 



+ 



-6s 



r = At + B 
9 



AA 2 



(3.20) 



where {A, B} are arbitrary constants. (We have perfomed a linear translation on t 
without loss of generality). It follows that the value of s is bounded from above by the 
constraint e 3s < 2v4/3 and we may therefore identify the eigenvalue of the momentum 
operator d/dr with the integration constant A, i.e. \p\ = 2A. 

In general, it is difficult to evaluate Eq. (|3.17|) exactly. However, it can be related 
directly to solution ( |3.12|) by performing a trivial rescaling c = —2c and specifying 



M(c) 



exp 



5 3 

h mc 

12 



(3.21) 



where m is an arbitrary constant. Substitution of variables ( ft.lOQ into Eq. fl3.17|) 
therefore implies that 



dc 



c c 



1/2 



exp 



-— + (u + v + 2m)- + — (a 
12 VP J 2 Ac 



(3.22) 



An integral of this form has been evaluated previously by Halliwell and Louko within 
the context of the path integral quantization of the four- dimensional de Sitter Universe 
41f] . They showed that it may be expressed in terms of products of Airy functions, 
where the specific combination is determined by the contour of integration. In our 
example, we wish to construct wavefunctions from a linear combination of bounded 
wavefunctions of the form (|3.2| ) and therefore require Rec < 0. If the contour of 
integration is chosen to lie along the line c = irj — e, where 77 is real and e > 0, it will 
pass to the left of the origin in the complex plane. In this case, Eq. ( |3.22j ) is given 

by 

* = Ai([M + m)Ai(u + m). (3.23) 
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Up to a numerical constant, this solution corresponds to Eq. fl3.12j ) with c\ = C3 and 
c 2 = c 4 = 0. It oscillates if either ^ + m<0or/i + m<0 and is exponentially 
damped when both arguments are positive. 

Thus far, we have derived exact solutions to the Wheeler-DeWitt equation. How- 
ever, it is useful to search for approximate solutions as well. Although such solutions 
are not exact, they can provide insight into the nature of the wavefunction. In the 
following Section, we shall discuss two classes of approximate solutions. 



4 Approximate wavefunctions 
4.1 Power series solutions 

Power series solutions to the Wheeler-DeWitt equation ( |3.1|) may be derived by ex- 
panding the wavefunction as the infinite sum of functions 



m=0 



This ansatz is a consistent solution to Eq. (|3.1|) if 



dad (3 







and 



d 2 * m 

dad (3 



vn — 1 ; 



m > 1. 



(4.1) 



(4.2) 



(4.3) 



Eq. ( |4.2|) is the canonical, one-dimensional wave equation and has the general solution 

% = P(a) + Q((3), (4.4) 

where P and Q are arbitrary, twice continuously differentiable functions. A modifi- 
cation of the Picard iteration scheme f2(| |34| may now be established by expressing 
\l/ m in terms of quadratures with respect to the null variables {a,j3}. When m = 1, 
Eq. (|4.3|) admits the separable solution 



(3 J da 1 P{a 1 ) + a j dfoQiPx 



(4.5) 



and this result may then be substituted back into Eq. (|4.3|) to derive \l/ 2 and so 
on. The general pattern is easy to deduce after a few iterations and we conclude, 
therefore, that 



* = ^0 + E 



-1} 



m=l 



ml 



+a r 



da r 



da-2 



daiPiai) 



03 



d(3 m ... d(3 2 d(3 l Q{(3 l ) 



(4.6) 
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is also a solution to the Wheeler-DeWitt Eq. (|3~T|). 

It should be emphasized that we have not assumed a semi-classical approximation 
in deriving these power series solutions. The advantage of this scheme is that the 
wavefunction is given in terms of arbitrary functions of a and ft. In many minisuper- 
space models these variables are related to the spatial volume of the Universe, where 
small values of a or ft typically correspond to small spatial volumes. Indeed, this is 
the case for two-dimensional dilaton-gravity cosmologies, since a is proportional to 
the square of the scale factor. Consequently, for P(a) = 0, we may view solution 
( (4.6| ) in the region of the origin as an expansion in powers of a small parameter a. 



4.2 Semi-classical wavefunctions 

Within the context of four-dimensional cosmologies, the nature of space-time is accu- 
rately described by classical physics when the spatial volume of the Universe is sig- 
nificantly larger than the Planck scale. It follows, therefore, that classical behaviour 
from the quantum regime should be predicted by the quantum theory. Presently, the 
problem of how such a transition might occur is an unresolved one. However, it is 
reasonable to suppose that the nature of semi-classical wavefunctions may provide 
some insight. 

In the WKB approximation, corresponding to the limit h — > 0, one treats some 
of the degrees of freedom {c} as classical variables and the remainder {q} quantum- 
mechanically. The wavefunction is then viewed as a linear superposition of waves 
of the form ^ « e - lS / h ^ where S is the classical action satisfying the Hamilton- 
Jacobi equation. This equation is derived by identifying the conjugate momenta in 
Eq. ( |2.18| ) with p a = OS /da and pp = dS/dft. It takes the form of a non-linear, 
first-order partial differential equation: 

dSdS 

where a new variable \ = v^ct has been introduced. It is well known that there 
exists a one-to-one correspondence between congruences of classical solutions and 



solutions to the Hamilton- Jacobi equation in two-dimensional minisuperspaces [f|2 
In principle, therefore, an arbitrary solution to Eq. ( |4.7|) may be generated once the 
classical solutions are known. 

However, parametric solutions may be found more directly by employing a Leg- 



endre transformation [Ell. We define new variables 



. dS OS 

and a new function 

P&V) = XZ + PV-S(X,P). (4.9) 
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Partial differentiation with respect to £ implies that 

x = |? = fr. ( 4 - 10 ) 

where the second equality follows from Eq. ( |4.7|) . Eq. ( |4.10| ) has the general solution 

p = \ve+f(v)i (4.H) 

where f{rj) is an arbitrary function of 77. 

A parametric solution may now be found by transforming back into the old vari- 
ables. After differentiation of Eq. (4.11) with respect to 77, we find that 



s = — fw + v-r- 

77 at] 



77^ arj 



It should be noted that this Legendre transformation is only self-consistent if the 
Jacobian 

c^Sc^S _ 

d^dfP ~ [d X df3 ) 

is non-vanishing. Solutions are said to be 'developable' if J 7^ and 'non-developable' 
if J = 0. All developable solutions can be written in the parametric form of Eq. (|4.12| ). 
In principle, we can determine rj = r](a,/3) from the second equation in ( |4.12| ) once 
the functional form of f{rf) has been specified. Substituting this result into the first 
equation yields the action in terms of the canonical variables, or equivalently, in terms 
of the original variables via Eqs. (|2.12|) , ( 2.15|) and ( 2.16 ). 



For example, the Jacobian is non-vanishing if 

(y ■ 

f(ri) = -+ftn, (4-14) 
77 

where /3j} are finite constants, and this ansatz leads to the action 

S = 2y/(a- ai ) (4.15) 



When / = 0, this solution is closely related to the exact solution of Eq. ( p.8|) that 
is given by * = (2y/a0), where (x) is the Hankel function. For small argu- 
ments this function has the asymptotic form 2iir~ 1 ln(x), so the wavefunction does not 
oscillate. On the other hand, the Hankel function has the form Hq\x) oc x~ l ^ 2 e ±lx 
at large arguments and this does have oscillatory behaviour. In this example, a large 
argument corresponds to a large value of the scale factor. Consequently, the argument 
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of this solution may be identified with the action (]4 .15 ) and represents a classically 
allowed solution that has tunneled from the Euclidean regime. 

Exact solutions can also be found if / oc r/ ±3 and / oc In 77. Furthermore, it is inter- 
esting to note that solution (|3.2| ) to the full Wheeler-DeWitt equation (|3.1| ) is of the 
WKB form \1/ = e~ tS , where S = ba + b~ x (3. This is an exact, non-developable solu- 
tion to the Hamilton- Jacobi equation (|4.7|) and in this sense the WKB approximation 
is exact for this solution. 

This concludes our discussion on approximate solutions. In the next Section we 
shall investigate whether other minisuperspace models can be solved in the manner 
discussed above. 



5 A class of integrable minisuperspaces 

It is interesting to investigate whether a wider class of models leads to the Wheeler- 
DeWitt equation Q3.1|). To proceed, we investigate an equation of the form 



d 2 



d 2 



dx 2 dy' 



+ 4m 2 (x, y) 



m = o, 



(s.i; 



where the superpotential, m 2 (x, y), is some function of the minisuperspace coordinates 
(x,y). 

We introduce new variables a = a(cr) and (3 = (3{r) that are arbitrary functions 
of the minisuperspace null coordinates a = x + y and r = x — y. These new variables 
satisfy the boundary conditions da/dx = da/dy and 0(3 /dx = —d/3/dy and these 
constraints ensure that the derivative terms in Eq. ( |5.1| ) are transformed into the 
canonical form: 

da df3 d 2 



dx dx dad/3 



+ m 



m = o. 



(5.2) 



It follows that Eq. (|5.2|) reduces to Eq. 
themselves solutions to the equation 



m 



da d(3 
dx dx 



if the new variables a and (3 are 



(5.3) 



da d(3 
da dr 



In principle, therefore, Eq. ( |5.1|) may be solved if a solution to the constraint equation 
Q5.3D can be found. Effectively, the problem of solving the linear, second-order partial 
differential equation Q5.1J ) has been reduced to finding a solution to the non-linear, 
first-order equation (|5.3|) and in many cases it is considerably easier to solve this 
latter equation. Indeed, it is clear from the second equality in Eq. ( |5.3| ) that when 
the superpotential has the generic form 



m 



m i 



(cr)m_ 



(5.4) 
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where m± are some known analytical functions, Eq. (|5.3|) admits the general separable 
solution a T 

xf da'm + (a'), (3 = A -1 / dr'm^r'), (5.5) 



a 

where A is an arbitrary separation constant. The region of minisuperspace covered 
by these coordinates is determined by the specific form of the superpotential. 

The Jacobian of the transformation leading to Eq. ( |5.2p vanishes whenever the 
null variables a = a(a) or f3 = f3{r) have turning points and these will occur at 
the zero points of the superpotential if Eq. (|5.3|) is satisfied. Thus, the Wheeler- 
DeWitt equation can be mapped onto the unit-mass Klein-Gordon equation if the 
superpotential is positive- or negative-definite over the whole region of minisuperspace 
covered by {a, (3} and, in addition, is a separable function of these null coordinates. 
For example, if a > and (3 < 0, the solutions discussed in Section 3.2 are also 
solutions to Eq. ( |5.1| ). If, on the other hand, the superpotential does vanish at some 
point in minisuperspace, equivalent transformations to those discussed above may be 
performed on both sides of the zero-point. The two solutions in the different regions 
may then be matched at the boundary by requiring that the wavefunction and its 
first derivative are continuous fl3| . 

There are a number of interesting minisuperspaces for which Eq. (|5.3|) can be 



solved exactly. In many cases the superpotential of the wavefunction is independent 
of one of the null coordinates, i.e. it is a single function of either a or r. This is the 
case for the Wheeler-DeWitt equation derived from a renormalizable, two-dimensional 
dilaton-gravity theory. One-loop quantum corrections to the CGHS action have been 
calculated by Russo, Susskind and Thorlacius [ 44| . In the conformal gauge g + - = 
9±± 



0, the one-loop effective action has the form 



S = - I d 2 x 



N 



-d +x d-x + -d + nd_n + ^- n ^ + W 9+fjd-fj 



(5.6) 



where 



X = W - - $ + e 



-2* 



represents a Liouville-type field, 



(5.7) 



(5i 



is a rescaled version of the dilaton field $, f) are conformal scalar fields and the 
constants k = (N — 24)/12 and /z 2 are assumed to be positive-definite^. 

2 The reader is referred to Ref. [Q for the details of the derivation. The numerical value of n is 
determined by including the one-loop contributions from the reparametrization ghosts, dilaton and 
conformal fields. The theory is one-loop finite if k = (N — 24)/12. 
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The Wheeler-DeWitt equation corresponding to this renormalizable model of 
dilaton-gravity has been derived by Mazzitelli and Russo []19 |. It has the form 



k d 2 k d 2 



N 



d 2 



4dQ 2 



^Edf 2 



4/Te 



2 2( Xo -!1o)/k 



^ = 0, 



(5.9) 



where xo ; etc., represent the zero modes of the harmonic expansion of the fields on 
the cylinder. In this analysis it is assumed that the coupling between the zero and 
higher-order modes is negligible and this is equivalent to invoking the minisuperspace 
approximation. This represents an improvement over the approximation employed to 
derive Eq. ( |3.1|) , however, since this latter equation follows from the classical action 
( f2.2| ), whereas Eq. ( [5.9D follows directly from the one-loop effective action ( [5.61 ). 

Eq. ( |5.9|) is solved by separating the wavefunction into its gravitational and 
matter components with the ansatz \& = &(xo, ^o)v 9 (/?o)- The plane waves <p = 



exp 



iJ2j Zjfjo 



form a basis for the solutions, where Zj are arbitrary constants. By 
identifying (x, y) = (x ,^o)> it is readily seen that $ statisfies Eq. ( |5.1| ), where the 
superpotential is given by 



Km 2 = ^_ ft ._ 2 _ 4 2 e 2( X o-fW« 

2 p 



(5.10) 



and Z 2 = J2j Zj represents the total momentum eigenvalue of the matter sector. If 
Z 2 < 2k + 4, m 2 is negative-definite and a function of (xo — ^o) only. We may 
therefore choose m + = 1 in Eq. ( |5.5[ ) and this implies that $ satisfies an equation of 
the form d 2 ^/dad(3 = — $, where 



1 



- « - 2 (xo 



fin 



a = A(xo + 
2 K (i 2 e 2(X0 ~ n ^ K 



On the other hand, the superpotential vanishes along the null line 

~Z 2 -2k-4 



Xo ~ 



In 



8/i 2 



(5.11) 



(5.12) 



if Z 2 > 2k + 4 and in this case different coordinate representations must be employed 
on either side of this line. 

We conclude, therefore, that the wavefunctions discussed in earlier sections also 
apply to this renormalizable model of dilaton-gravity. In particular, Lorentzian solu- 
tions to Eq. ( |5.9| ) may be generated from linear superpositions of Euclidean solutions 



and vice-versa when a > and (3 < 0. It follows immediately from Eq. (|5.11|) that 
these conditions are satisfied for all xo > > when A > and Z 2 < 2k + 4. 

This concludes our discussion on exact bosonic wavefunctions. In the following 
Section we shall investigate whether supersymmetric extentions to the quantum mod- 
els discussed above can be performed. 
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6 Supersymmetric quantum cosmology 



Graham discovered that a hidden symmetry exists in the Bianchi IX Universe by 
showing how the classical superspace Hamiltonian may be viewed as the bosonic 



part of a supersymmetric Hamiltonian [45|. It has now been shown that this hidden 



symmetry exists in all Bianchi class A models [|!| |4^], ^8], This implies that a 
supersymmetry can be introduced at the quantum level. This supersymmetric exten- 
sion of the quantum theory has significant consequences for quantum cosmology, as 



shown by calculations for the Bianchi II Universe [|9|. It is thought that these exten- 
sions may provide valuable insight into some of the questions relevant to a complete 
theory of quantum gravity. In particular, they may resolve the problems encoun- 
tered when one attempts to construct a conserved probability from the wavefunction 
p8| . It is therefore of interest to investigate whether hidden supersymmetries exist 
in two-dimensional dilaton-gravity models. 



We begin by briefly reviewing the 'hidden symmetry' method of Graham H5" 



In the minisuperspace approximation the classical Hamiltonian constraint takes the 
form 

2H = G^ptfv + W(q) = 0, (6.1) 

where G^ u is the metric with signature (— , +,+,.. .) on the (D + l)-dimensional min- 
isuperspace spanned by the finite number of degrees of freedom q^ (// = 0, 1, . . . , D). 
The momenta conjugate to these variables are p M and W represents the superpoten- 
tial. This Hamiltonian is the bosonic component of a supersymmetric Hamiltonian 
51] if there exists a function I(q) that respects the same symmetries as H and 



is itself a solution to the Euclidean Hamilton- Jacobi equation 

dl dl 

W = (6.2) 
dq* dq u v ' 

Fermionic degrees of freedom ip^, (p v obeying the spinor algebra 

[yA y/] + = r , ^] + = 0, ^] + = CT (6.3) 

are then introduced. It follows that the supercharges 

satisfy 

Q 2 = Q 2 = (6.5) 
and, if Eq. ( |6.2| ) is satisfied, the Hamiltonian ( |6.1| ) may be written as 

2H = [Q,Q] + , [H Q ,Q\- = [H ,Q]- = 0. (6.6) 

These equations represent the algebra for a single, complex supersymmetry charge Q 
and the model therefore exhibits an iV = 2 supersymmetry |)T] . 
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This symmetry is preserved at the quantum level by choosing the representation 
(pP = Qt 1 and ^ = G fJ-X d/d6 x for the fermionic degrees of freedom, where 9^ are Grass- 
mann variables [[S [IT] . The bosonic degrees of freedom have the usual representation 

= —ihd/dq^ 1 . Eqs. (|6.5|) and ( |6.6|) now represent the operator realizations of the 
super symmetric algebra. The quantized superspace Hamiltonian is given by 

H = H ^\wk^' v " ] - (6 ' 7) 

and has an additional term that vanishes in the classical limit. The existence of 
this term suggests that suitable imaginary or complex solutions to Eq. (|6.2|) will be 
difficult to find. It follows that the supersymmetric wavefunctions are annihilated by 
the supercharge operators 

Qty = = (6.8) 

and it is these constraints that represent the 'square roots' of the Wheeler-DeWitt 
equation. 



6.1 Induced gravity theory 

To investigate whether the two-dimensional cosmological models considered in Section 
2 exhibit a hidden supersymmetry of the form discussed above, we must first identify 
the symmetries of the classical Hamiltonian (|2.14 ). The kinetic part is invariant 



under the simultaneous interchanges u <-» ±t> . However, the full Hamiltonian is not 
necessarily invariant under this interchange because of the generality of the dilaton 
potential. On the other hand, it is symmetric under both simultaneous interchanges 
if the potential is constant, i.e. if V = A 2 . This form of the potential in Eq. ( |2.8j ) 
corresponds to the induced gravity action for q = and it is straightforward to 
verify that the above symmetries are equivalent to an invariance under a «-> ±/3. 
Evaluation of Eq. ( [2.16|) implies that (3 = q 2 \ 2 e~^^ q and comparing Eq. ( |2.18|) 



with Eq. (|6.1|) implies that the non-zero components of the minisuperspace metric 
are G a p = Gp a = a\ 2 . The superpotential is therefore given by W = —2a\ 2 and it 
follows that one solution to Eq. (|6.2| ) that respects the symmetries of the Hamiltonian 
is 

/ = -2i(af3) 1/2 . (6.9) 

Since we require this 'Euclidean' action to be real, we must assume that A 2 < 0. 

The functional form of the supersymmetric wavefunction may now be determined 
by solving the constraints ( |6.8| ). Due to the anticommutation relations obeyed by the 
Grassmann variables 8^, the general supersymmetric wavefunction may be expanded 
as 

^ = A+ + B 6° + B x 9 l + C 2 9°e\ (6.10) 

where the bosonic functions {A + , B , Bi, C 2 } are functions of {a, 0} only. The anni- 
hilation of the wavefunction by the supercharge operators then translates into a set 
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of coupled, first-order partial differential equations: 







d + dl' 

da da 










' d df 

_dj3 + d(3_ 






a or 

da da 




' d dl 

_d(3 + d(3 


d dl' 


B,+ 


' d dl 


da da 


_df3~d[3 
' d dl 
da da 
d dl' 

dj3~d(3 



A, 

A, 








B = 



B 
C 2 
C 2 ~- 



-- 
= 
0. 



The solution to these equations is given by 

A+ = e- 1 

_ dF p dl 

da da 
_dF dl 



where the function F = F(a,/3) is itself a solution to the equation 

F = 0. 



d 2 F ( d 2 I dl dl\ 

dad (3 \ dad (3 da d/3 J 



When / is given by Eq. Q6.9Q , Eq. (|6.13|) simplifies to 

d 2 d 2 



dw 2 dz 2 



w z + z z + 2 



0. 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



where {w,z} are defined by Eq. ( |2.2(J| ). Hence, F may be interpreted physically as 
the wavefunction for a quantum system describing two coupled harmonic oscillators 
that have identical frequencies but a difference in energy of 2. The solution to Eq. 
( |6.14j ) has the separable form 



F = H n (^ftf + H n+l - e 



3 -(cH-7/3) 



(6.15) 



The functions A + and C 2 represent the empty and filled fermion sectors of the 
Hilbert space. Both may be interpreted as lowest-order, WKB approximations to 
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exact solutions of the bosonic Wheeler-DeWitt equation ( |3.8| ). This equation may be 
written formally as H(o)^ = Hm^f, where we have split the Wheeler-DeWitt operator 
into the two components 

H(o) = ~ ~ 97e 6s , H W = ~ (6.16) 

Application of Hm implies that H^ty = E^f, where E = p 2 . If the eigenvalue 
p is real, this equation may be interpreted as the Schrodinger equation, where E 
represents the energy associated with ||52|| . Therefore, the state with E = 



corresponds to the state of lowest energy. When (3 < and p — 0, the general form 
of the bosonic wavefunction (|3.9| ) is given by a linear combination of modified Bessel 
functions Iq(x) and Kq(x), where x = 2^a^f3. For large x these functions have the 
asymptotic forms Jo oc e x and K oc e~ x , respectively, and these limits correspond 
to the solutions C<i and A + . The supersymmetric vacua are therefore closely related 
to their semi-classical limits and correspond to the pure bosonic states of lowest 
energy. These features appear to be generic properties of supersymmetric ground 
state wavefunctions ET . 



6.2 A string inspired model 



String theory exhibits a symmetry known as target space duality [[53]. (For a recent 
review see, e.g., Ref. [0). In two-dimensional space-times, a string cannot tell if 
it is propagating on a circle of radius a or of radius a~ x . In effect, this allows one 
to transform between theories of radii a and a~ l after a suitable translation on the 
dilaton field j55|. It is convenient to consider this symmetry within the context of 
the action fl2.2p in the absence of loop corrections, i.e. D = c > 0. The Hamiltonian 
derived from this action for the cosmological space-time ( |2.10| ) is given by 



H 



-2$ 



a®' 



CM 



(6.17) 



and it is straightforward to verify that Eq. ( 6 . 1 7| ) is invariant under the duality 
transformation 

a = -, $ = + lna. (6.18) 



If we introduce the coordinate pair 

X = ae"*, 
the Hamiltonian takes the form 



Y 



(6.19) 



H = -—XY-cXY 



(6.20) 
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and invariance under the duality transformation Q6.18 ) is therefore equivalent to an 
invariance under the simultaneous interchange X Y . 

The momenta conjugate to X and Y are px = 4Y/N and py = AX/N, re- 
spectively. It is convenient to perform a rescaling of these degrees of freedom by 
defining a = X 2 and (3 = Y 2 . The classical Hamiltonian ( |6.20| ) is then given by 
Eq. ( |6.1|) , where the non- vanishing components of the minisuperspace metric are 
G a p — Gp a = (a/?) 1 / 2 and the superpotential W = 2c{a(3) 1 ^ 2 . There exists a hidden 
supersymmetry if / satisfies 

dl dl . « 

T a w = c (6 ' 21) 

and also respects the duality symmetry a «-> /3. One solution satisfying the necessary 
conditions is / = 2y/ca(3 and since {a, j3} are positive-definite, c > is necessary for 
the solution to be real. 

We conclude, therefore, that the general supersymmetric wavefunction may also 
be found in closed form for this theory. The bosonic functions in the Grassmann basis 
expansion of the wavefunction are again given by Eq. ( |6.12j ), but F has the slightly 
different form 

F = H n ( V )H n+1 (Oe-^ + ^ /2 , (6.22) 
where £ = c 1 ^^ 1 ' 2 + (3 1 / 2 ) and rj = c x l\a x l 2 - /3 1 / 2 ). 



7 Conclusions and discussion 

In this paper we have investigated the quantum cosmology of a generalized class of 
two-dimensional dilaton-gravity models. If the dilaton potential contains no roots, 
i.e. if V(ip) 7^ for all physically interesting x/j, the classical dynamics of these 
Universes is equivalent to that of a non-interacting, point particle propagating over a 
portion of two-dimensional Minkowski space. A large subset of this class of models is 
dynamically equivalent to the isotropic, constrained oscillator- ghost-oscillator system. 
This suggests that the relationship between quantum configurations and classical 
space-times, as discussed in Ref. |17| , could be generalized to these models. 

Furthermore, these correspondences imply that the Wheeler-DeWitt equation can 
be expressed as the unit-mass Klein-Gordon equation if a suitable choice of factor 
ordering is made. This allows a number of exact and approximate solutions to be 
found. Quantum states corresponding to Lorentzian geometries may be generated 
from an infinite sum of Euclidean solutions and vice-versa. The Hamilton- Jacobi 
equation can be solved by employing a Legendre transformation and all developable 
solutions to this equation were found in parametric form. 

We proceeded to identify a wider class of integrable two-dimensional minisuper- 
spaces that can be solved exactly by mapping the Wheeler-DeWitt equation onto the 
unit-mass Klein-Gordon equation. This mapping is possible if the superpotential of 
the wavefunction is a separable function of the null coordinates over minisuperspace. 
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We applied this result to the Wheeler-DeWitt equation derived from a renormalizable, 



two-dimensional dilaton-gravity model [19 



One of the main problems with the quantum cosmology program is the construc- 
tion of a non-negative norm from solutions to the Wheeler-DeWitt equation. This 
equation is a hyperbolic, second-order partial differential equation, so the conserved 
current associated with it is not necessarily semi-positive definite. Consequently, it is 
not clear that such a current will provide a suitable measure of probability. A similar 
problem is encounted when the Klein-Gordon equation for a scalar field is solved. In 
this case, however, the ambiguity is resolved by taking the 'square root' and in view of 
the close correspondence between the Wheeler-DeWitt and Klein-Gordon equations, 
it has been suggested that a similar technique might solve the corresponding problem 
in quantum cosmology [[f^, [5(| . 

This suggests that one should search for supersymmetric extensions to quantum 
cosmology. It was shown in Section 6 that the classical Hamiltonians derived from 
the induced gravity theory and a string-inspired model may be viewed as the bosonic 
components of a supersymmetric Hamiltonian. In the latter case, the origin of this 
symmetry may be traced to the invariance of string theory under duality transfor- 
mations. The hidden symmetry method was employed to derive the corresponding 
quantum constraints for the two models. This method differs from other approaches 
to supersymmetric quantum cosmology because it does not start from a field theory 



of supergravity [ |j| . The quantum constraints can be solved exactly and closed-form 
expressions for the general supersymmetric wavefunction were found. It would be in- 
teresting to investigate whether this method can be applied to more general models. 

An alternative approach to quantum cosmology is the third quantization proce- 
dure ^7\, The aim of this approach is to construct a consistent probabilistic 
measure in quantum gravity by promoting the wavefunction of the Universe to a 
quantum field operator that acts on a Hilbert space of states. The 'vacuum' state 
in this space is identified as the state where the Universe does not exist. Topology 
changing processes can then be described by including self-interactions of the Uni- 
verse field. Moreover, in the minisuperspace approximation a suitable combination 
of the dynamical degrees of freedom may be associated with a time variable in the 
Wheeler-DeWitt equation. It then follows that the superpotential of the wavefunc- 
tion may be viewed as a 'time-dependent' function. In ordinary quantum field theory 
it is well known that particles are created from the vacuum by a time-varying exter- 
nal potential and this suggests that Universes could be created via a similar process. 
In practice the Universe field is expanded into positive frequency in- and out-mode 
functions and their hermitian conjugates. The in- and out-modes are related to one 
another by the Bogoliubov coefficients and these determine the number of Universes 
in a given mode p9| . The creation of Universes in this picture arises because the two 
Hilbert spaces generated by the in- and out-mode functions are inequivalent and this 
results in non-zero Bogoliubov coefficients. 

Recently Vilenkin |60|] has argued against this picture of Universe creation. His 
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main objection is that the time variable constructed in minisuperspace models is 
generally not a monotonically increasing function since Universes can contract as 
well as expand. He then interprets the creation of a pair of Universes in terms of a 
contracting Universe that begins reexpanding at a finite radius. On the other hand, 
he does suggest that third quantization might be appropriate for describing topology 
changing processes in two-dimensional Universes. 

Since there is currently no generally accepted interpretation of third quantization, 
it is of interest to investigate its consequences further. The procedure can be applied 
to the class of two-dimensional models ( |2.£| ) for which (3, as defined in Eq. ( |2.16| ), 
is positive-definite for all ip. The variables (|3.7| ) take all values in the range (s, r) G 
(—00, +00) and we may therefore view s as the time variable in the Wheeler-DeWitt 
equation ( |3.8| ). The scale factor of the Universe vanishes as s — ► — 00 and infinite 
spatial volume corresponds to the limit s — > 00. Formally, this model is identical to 
the one considered previously by Hosoya and Morikawa |58| , so their results will apply 



here. The appropriately normalized positive-frequency in- and out-mode functions are 
given by 



r) = I' V sinh^f e*»7„(e-) 




1/2 



e -n\p\/6 e ipr H (2) / 3^ ? ( ? ^ 



respectively, where v = —i\p\/3. As s — > oo, u° nt oc e tS and this is the WKB solution 
given by Eq. ( 4.15| ) with / = 0. It follows that the Bogoliubov coefficients are given 



by ci{p, q) = b i$pq and c 2 (p, q) = b 2 5 pq , where 



bl = (l- e- 2 ^)~ 1/2 , \b 2 \ = (e 27 ^/ 3 - lY 1 ' 2 (7.2) 



and the average number of Universes with 'energy' p therefore has a Planckian dis- 
tribution 

iV P =|c 2 (p,p)| 2 =(e 2 ^ l/3 -l) _1 . (7.3) 

Hosoya and Morikawa |5£| extended this free field theory by including a \1/ 3 in- 
teraction that describes the spliting of a 'mother' Universe into two 'baby' Universes 
of identical topology. By treating the mother Universe in a classical fashion, they 
showed that the quantized baby Universes also have a Planckian distribution. The 
formal equivalence of their model with those studied in this work suggests that similar 
conclusions should apply for a wide class of two-dimensional cosmologies. It would 
be of interest to investigate these possibilities further. 
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